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ABSTRACT 

Progenitors of Type la supernovae (SNe) have been predicted to modify their ambient circumstellar 
(CSM) and interstellar environments through the action of their powerful winds. While there is X-ray 
and optical evidence for circumstellar interaction in several remnants of Type la SNe, widespread 
evidence for such interaction in Type la SNe themselves has been lacking. We consider prospects for 
the detection of CSM shells that have been predicted to be common around Type la SNe. Such shells 
are most easily detected in Na I absorption lines. Variable (declining) absorption is expected to occur 
soon after the explosion, primarily during the SN rise time, for shells located within ^ 1-10 pc of a 
SN. The distance of the shell from the SN can be determined by measuring the time scale for line 
variability. 

Subject headings: ISM: bubbles — supernovae: general 
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1. CSM SHELLS AROUND TYPE lA PROGENITORS 

Type la supernovae (SNe) play an important role in 
many areas of astrophysics, even though their progenitor 
systems and detailed explosion mechanisms are poorly 
understood. One fundamental problem involves the pre- 
diction from many models, but so far nondetection, of 
substantial circumstellar material (CSM), likely in the 
form of shells at pc-scale distances. Here, we predict the 
decrease in Na I absorption at early times that such shells 
should produce, a companion effect to Na absorp tion in- 
creases seen somewhat later (jPatat et al.ll2007a l). 

Most popular current Type la scenarios involve sin- 
gle degenerate ("SD") progenitors, where the donor star 
supplies mass to the white dwarf (WD), driving it close 
to the Chandrasekhar limit, where instabilities produce 
an explosion. A crucial step in this process is the ability 
of an accreting white dwarf to expel material accumu- 
lating on its surface in excess of a critical accretion rate 
through a fast {vw ~ 1000 km s~^) radi atively driven 
stellar outflow (|Hachisu et al.lll996l . fl999afH ). In the ab- 
sence of such a wind, formation of a common envelope 
around the progenitor system is unavoidable because of 
the expansion of the WD envelope. The companion then 
merges with the white dwarf, preventing an explosion. 

The WD accretion winds are expected to have pro- 
found consequences on the CSM and interstellar medium 
(ISM) around Type la SNe. Fast and powerful winds 
not only sweep denser CSM out of the progenitor sys- 
tem, they have also been suggested to strip m aterial 
from the donor star (e.g., iHachisu et al.l HQQI) which 
then could form an asymmetric disk in the binary orbital 
plane. When swept by the SN blast wave, this material 
would be expected to produce high velocity absorption 
features, and perhaps detectable radio and X-ray emis- 
sion. Farther out, the white dwar f wind blows up large- 
scale bubbles in the ambient ISM. iBadenes et al.l ()2007f ) 
investigated the evolution of these bubbles with one- 
dimensional hydrodynamical simulations. For an ISM 
hydrogen density no of 0.43 cm^"^ (po = lO"^'* g cm^'^) 
and Vyj = 10'^ km s^^, large (> 15 pc) bubbles are blown 
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by progenitor winds, sweeping the ambient ISM into mas- 
sive (> 300 Mq) shells. 

Since typical progenitor space velocities of 20 
km s~^ are larger than final shell velocities, a SN ex- 
plosion generally occurs after the progenitor has moved 
away from the bubble's center, by a distance depending 
on no and on the progenitor scenario under considera- 
tion. If the WD explodes early ('^ 5 x 10^ yr) after 
th e onset of the WD w ind (as in models HI and HP3 
of IBadenes et al.ll2007D . the explosion occurs off-center 
but still inside the bubble. For no = 0.43 cm~'^ and 
= 20 km s~^, the bubble's radius is 17 (22) pc for 
model HI (HP3), while the SN exploded 7 (9) pc away 
from the shell. The SN would be deeper inside the bub- 
ble for slower moving progenitors and less dense ambient 
ISM. For faster progenitors and denser ISM, an explo- 
sion could even occur outside the bubble, but still in the 
vicinity of the swept-up shell. 

Stellar motions are most important for progenitors that 
explode long after the onset and cessation of the wind 
(2 X 10^ and 10^ yr, re spectively, in models HPl and L2 
of IBadenes et al.l 120071 ). A progenitor with = 20 km 
s""'^ travels 40 pc in 2 x 10^ yr, beyond bubbles with a 
radius of 27 (36) pc in model HPl (L2) (for no = 0.43 
cm~^). Without a fast wind blowing at the time of the 
explosion, these SNe explode in the undisturbed ISM. 

Another outcome is expected for a progenitor with 
a moderate, long-duration wind that is still blowing 
when the SN explo des (as in models LI and LVl of 
IBadenes et al.l l2007l ). Even if the progenitor left the 
bubble's interior, the wind interaction with the ambi- 
ent ISM leads to the formation of a bow shock ahead 
of the moving SN progenitor (e.g.. Iwilkinl Il996[ ). The 
standoff distance Rq of the swept-up ISM shell in the di- 
rection of the stellar motion is 0.75 pc for a mass-loss rate 
= IO-'^Mq yr-i, - 10^ km s"!, = 20 km s-\ 

and no = 1 cm~'^ (Rq cx Mw'^vU^nQ^^^v^^). The dis- 
tance to the shell is 1.3 pc in the plane perpendicular to 
the stellar motion direction and passing through the pro- 
genitor. The shell surface density is ~ Rono (= 2.3 x 10^^ 
cm~^ for Rq = 0.75 pc and no = 1 cm"^; I WilkinI [l996l : 
IComeron fc Ka pcr 1998). As discussed next, such a CSM 
shell and more massive CSM shells swept up by the 
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WD winds are expected to produce detectable absorp- 
tion lines in SN spectra. 

CSM has also been revealed by young Type la super- 
nova remnants (SNRs). Kep ler's SN was a Type la ex- 
plosion within a dense CSM (iRevnolds et al.ll2007l ). The 
total shocked CSM mass in Kepler's SNR is ~ 0.75Mq 
air et al.ll2007[ ). The blast wave propagated to 2 pc 
away from the explosion center, and it is still encounter- 
ing dense material with a preshock density of ~ 5 cm~^. 
This CSM was ejected by the SN progenitor, unlike the 
massive swept-up ISM shells just discussed. We generally 
expect shells at small radii to be dominated by ejected 
CSM, with the swept-up ISM dominating at large radii. 

2. A SIMPLE CSM SHELL MODEL 

A geometrically thin, spherically symmetric CSM shell 
with radius Rg and mass Ms has a hydrogen column 
density 
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N{H) = 1.78 X 10^ 



Rs 
2 pc 



For a shell swept up by a progenitor WD wind in a uni- 
form ISM with density no, equation ([1]) becomes 



N{H) = 2.06 X 10 



18 



Rs 



no 



(2) 



2 pc cm ^ 

Theoretical considerations and observational evidence 
suggest a wide range in shell masses, from a fraction of 
a solar mass at sub-pc scales to hundreds of solar masses 
at ^ 10 pc scales (Section [1]). This simple CSM shell 
model does not encompass progenitors that left the am- 
bient ISM undisturbed. 

Shells with N(H) - lO^^-lO^^ cm'^ (eqs. fflandE]) can 
be most easily detected in optical Na I and Ca II ab- 
sorption lines. We focus on Na I lines, because of the 
likely depletion of Ca onto dust that may reduce Ca II 
absorptio n below detectabl e levels. The CSM dust found 
m Kepler HI air et al.ll2007[ ). for example, would result in 
significant depletion of Ca. Neutral Na is quickly ionized 
by the Galactic UV radiation field because of its low (5.14 
eV) i onization potential. For the "standard" UV field 
(e.g.. lPequignot fc Aldrovandilll986f ). the ionization rate 
r(Na I) is equal to 1.0 x 10~" s~^ As in the general ISM, 
the Na ionization state is determined by the balance be- 
tween ionizations and recombinations, n(Na I)r(Na I) = 
nen(Na II) aNaii{Te), where a^aii is the recombination 
coefficient for Na II. For gas with predominantly singly 
ionized Na with a solar (cosmic ) abundance of 2.0 x 10^^ 
(by number with respect to H. iLodders et al.l 2009) and 
"standard" Galactic UV radiation field, N(Na I)/N(H)== 



2.2 X 10"^ieaNaii(Te)/aNaii(1037^), and 
7V(NaI) = 

I QNall(T'e) He / Rs 

a^^ii{10^K) 0.1 cm-3 ^2 pc 



3.9 X 10' 



Mr, 



For a shell swept up by the progenitor WD wind, 



7V(NaI) 
4.6 X 10 



(3) 



(4) 



10 aNall(re) Ue 



Rs no 



aNaii(lO'^-ftr) 0.1 cm~3 2 pc cm^^ 

The electron density (and to a lesser degree also the 
electron temperature Tg) within the CSM shell are ad- 
ditional parameters that control the abundances of trace 




Fig. 1. — Normalized Na I ionization rate as a function of time 
since explosion (solid line: based on spectral templates of Nugent 
et al. 2002; dotted line: on the SNLS templates). A synthetic Swift 
uvw2 light curve (dashed line; based on Nugent et al. 2002), and 
normalized uvw2 fluxes of SN 2007af (stars) are also plotted. 

ions such as Na I and Ca II. Here, some degree of ion- 
ization of H produces the free electrons that permit the 
existence of neutral Na. A modest photon ionizing flux, 
provided by a SN progenitor or by the ambient extreme 
UV Galactic radiation field, is likely to ionize a small 
(r-^ 0.1) fraction of H atoms, so the compressed shell gas 
at typical densities of ~ 1 cm^'^ (Badcnes ct al. 2007() is 
expected to be detectable in Na I lines. 

In the SD progenitor scenario, ionizing photons are 
produced by t he WD and its accretion disk. For example, 
IChugail ()2008f ) considered a blackbody with T = 3 x lO" 
K and L = 300^© as an ionizing source in his con- 
siderations of absorption lines in SNe with symbiotic 
progenitors. Such a modest luminosity source is ca- 
pable of significant ionization of the ambient gas to 
large distances away from the progenitor. We demon- 
strate this by performin g photoionizat ion calculations 
with Cloudy (ver. 08.00: fl^rland et al.l ll9981. In addi- 
tion to the blackbody radiation produced by the progen- 
itor, we also include the diffu se Galactic UV radiation 
field (Drai ne fc Bertoldilll996 !) that is responsible for the 
ionization of Na I. A CSM sheU with 1 Mq, i?^ = 2 pc, 
and n = 5 cm~^ becomes nearly completely (99%) ion- 
ized by the progenitor radiation. But Na remains singly 
ionized, leading to an appreciable (5 x 10^° cm~^) Na I 
column density. For a more massive (10 Mq) and dis- 
tant (10 pc) shell, the H ionization decreases to 80%; 
N(Na I) drops slightly to 4 x 10^° cm-^. N(Na I) in- 
creases with the increasing shell mass; for 50 Mq, it is 
equal to 2 x 10^^ cm~^, while H becomes partially (49%) 
neutral. There are large spatial variations of the diffuse 
ionizing radiation in galaxies, and its penetration into 
the CSM shell depends strongly on the assumed geome- 
try and dynamical history of the shell. In view of these 
large uncertainties, we simply use equations (3) and (4) 
in our estimates of Na I column densities, using rie as a 
free parameter to describe the unknown ionization frac- 
tion of H in the CSM sheU with density ~ 1-10 cm"^. 

3. VARIABLE NA I ABSORPTION 

Photoionization of Na I by SN photons is very effi- 
cient in th e immediate (< several xlO^^ cm) vicinity 
of the SN (Chugai 2008), decreasing rapidly with dis- 
tance. We considered photoionization of Na I as a func- 
tion of shell location Rs and time t since the explo- 
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Fig. 2. — Decrease in Na I absorption as a function of time 
and distance from the SN. The curves are labeled by the amount 
of decrease in N(Na I). Reasonable prospects for the detection of 
absorption line variability occur within the cross-hatched area. 

sion, using Type l a Branch-normal spectral templates 
(jNugent et al.ll2002D^ and Na I photoionization cross sec- 
tions from IVerner et all (jl996l ). Only photons with en- 
ergies larger than the Na I ionization potential of 5.14 
eV (A < 2410A) are capable of ionizing neutral sodium. 
Figure [1] shows the resulting Na I ionization rate as a 
function of time, normalized to 1 at its maximum. Wit h 
My = -19.46 at maximum hght (jTammann et al.ll2008D . 
the normalizing factor is equal to 3.3 x 10~^ (i?s/2 pc)~^ 
s~^. We also show in Figure [1] the normalized Na I ion- 
ization rate based on spectral templates from t he Su- 
pernova Legacy Survey fSNLS: [Hsrao et al.|[2007l )'^. The 
normalizing factor is only 7% smaller than for Nugent 's 
templates. Flux measurements with the Swift satellite's 
uvw2 filter [Xcentrai = 1928A, FWHM= 657A) are of 
particular interest, since they provide light curves short- 
ward of the Na I ionization limit. Aside from a notice- 
able excess at late times caused by a red leak in this 
filter, a synthetic uvw2 light curve matches the Na I 
photoionization rate well (Figure [T]). At its maximum 
on day 15, M,„„„2 = -17?' 5 (using Swift UVOT cali- 
bration bv iPoole et al.l 120081 ) . For comparison, we also 
plot a normalized uvw2 light curve for SN 20G7af (from 
iBrown et al.ll2009f l. a nearby, well-observed normal Type 
la S N. We assumed a typical Type la rising time of 18 
days (|Garg et al.ll2007l) in the B band for SN 2007af. The 
calculated Na I ionization rate curve is narrower than 
the uvw2 light curve of SN 2007af by about 3 days, while 
the total number of Na I ionizing photons is likely to 
be reasonably well modeled by the template (based on 
a dereddened Muvw2 of — 16™8 for this SN). This mod- 
est mismatch likely reflects systematic deficiencies in our 
current knowledge of Type la UV spectra, due to their 
UV faintness, uncertain reddening corrections, and in- 
trinsic UV spectral variations. 

Figure [2] shows how Na I absorption decreases as a 
function of Rs and t, using the spectral templates of 
iNugent et al.l (|2002f ). Because the Na I ionization rate 
peaks just before the i?-band maximum at ~ 18 days 
(Figured]), variable absorption is expected in the first 
20-30 days after the explosion, primarily during the SN 
rise time. Detection will be difficult for decreases below 

^ Available online at www.supernova.lbl.gov/~nugent 
^ Available online at www.astro.uvic.ca/~hsiao 
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Fig. 3. — Remaining Na I fraction as a function of time since 
explosion for shells located at 1, 2, 5, and 10 pc (solid lines from 
left to right) from the SN. The dashed curve corresponding to the 
uniform ambient ISM is also shown. 

10% or where 99% of Na ions are ionized. Prospects for 
detection of variability are optimal in the middle of the 
cross-hatched area in Figure [21 and observational con- 
straints favor detection of shells at pc-scale distances. 

Figure [3] shows the remaining Na I fraction for shells 
with Rs = 1, 2, 5, and 10 pc. These curves depict vari- 
ations in N(Na I) (or the equivalent width of the Na D 
lines in the limit of low line optical depths) as a function 
of time. The rate of change in absorption is much faster 
for shells at smaller Rs, reflecting distance-dependent 
(oc i?7^) variations in the Na ionization rate. This sen- 
sitivity to distance from the SN can be utilized to infer 
distances of CSM shells from observations of absorption 
line variability. We also show in Figure [3] the remaining 
Na I fraction for uniformly (up to 10 pc) distributed ISM. 
Except at early times, variations in absorption are rather 
similar to variations produced by the shell with Rs = 5 
pc. But absorption lines from an undisturbed ISM will 
be harder to detect than from the CSM shells, because 
these lines are more likely to be located within saturated 
cores of the Na I D absorption line profiles produced by 
the ISM in host galaxies. 

Intrinsic variations in the UV among SNe are likely 
to affect the detectability of Na I lines. Overluminous 
SNe such as SN 1991T are attractive targets, because 
their higher than average fluxes facilitate high spectral 
resolution observations of absorption lines, while their 
potentially higher UV luminosities may allow probing 
of the swept-up ISM to larger distances. The bright- 
est S Ne are also associa ted with star-forming galaxies 
fe.g.. lSullivan et al.|[2006l ). implying more gas-rich envi- 
ronments near overluminous SNe. Conversely, sublumi- 
nous SNe such as SN 1991bg are more difficult to observe, 
and their potentially lower UV luminosities might allow 
probing of the CSM only to sub-pc distances from the 
SN. 

4. DISCUSSION 

Variable Na I absorption in Type la SN spectra probes 
CSM at pc-scale distances from the SN, such as the CSM 
seen in Kepler's SNR. With i?, = 2 pc, Ms = 0.75Mq, 
and Ue — 0.5 cm~'^, the preexplosion N(Na I) is 1.5 x 10^^ 
cm^^ (from eq. 3 with — 10"^ K). The decrease in ab- 
sorption is shown in Figure [3] (second curve from left); 
N(Na I) drops below 10^° cm^^ ^g^y 20. Sufficiently 
rapid motion (radial velocity of —230 km s~^: lBlair et al.l 
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Il991h of Kepler's SN progenitor would have assured no 
overlap with the (presumably) saturated ISM absorption 
in the line of sight to Kepler's SN. It is possible to mea- 
sure the predic ted Na I D absorpti on in Type la spectra 
at early times. iPatat et"an (|2007bf ) obtaine d a stringent 
2 X 1 0^° cm-2) upper limit for SN2000cx. ISimon et all 
2007f) observed SN2007af at three epochs (days 14, 35, 
and 42 in Figure [T]). Their observation on day 14 was 
the least restrictive; the Sct limit on the CSM absorp- 
tion is 9 X 10^° cm~^. The predicted column density on 
day 14 is 5 X 10^" cm^^, at the detection limits of this 
observation. No detection would have been expected on 
days 35 and 42. Our simple shell model for Kepler is 
of course uncertain. Kepler's SNR is very asymmetric; 
depending on viewing angle, much lower or much higher 
column densities are possible in a SN with such strongly 
asymmetric CSM. Kepler's SN also has a central emis- 
sion enhancement, indicating the presence of CSM at 
distances less than 2 pc from the SN; that would result 
in faster ionization by the SN. 

The remnant of SN 185, SNR ROW 86, prob- 
ably res ulted from a Type la explosion within a 
bubble (iBadenes et all 120071 ). At -2.5 kpc away 
(jSoUerman et al.ll2003f) . its 21' radius corresponds to 15 
pc. Assuming a spherically symmetric bubble with this 
radius, the mass of the swept-up ISM is equal to 5G07io 
Mq. From Equation (4), the preexplosion N(Na I) is 
larger than 10^^ cm^^ for uq > 0.29 cm~^ (with rie = 0.1 
cm"-^ and Tg — 10'^ K). RCW 86 is much brighter in the 
SW than elsewhere, suggesting an off-center explosion in 
this region of the SNR. The offset of the SN progeni- 
tor from the bubble's center is not known; we assume 
— 10 pc offset from the center (~ 5 pc distance from the 
shell) . The decrease in absorption is shown in Figure [3] 
(second curve from top); half the Na atoms are ionized 
in the shell section closest to the SN. Significant (> 10%) 
variations in N(Na I) would be seen by about half of ran- 
domly distributed observers; for the other half, the CSM 
shell would have been too far away from the SN to pro- 
duce variable absorption. The predicted column densities 
depend sensitively on geometric details, but our simple 
model for RCW 86 demonstrates that shells swept by SN 
progenitors should be detectable in Na I lines. 

The variable Na I absorption d etected in SN2006X 
(|Patat et al.ll2007a[ ) and SN 1999cl fBlon din et 311120091) 
cannot be due to our CSM shell effect, as it is stronger, 
slower, and increasing rather than decreasing. While the 
origin of these variations is still not understood, both 
SNe are among the niost hig hly reddened Type la SNe, 
leading [Blondin et al.l (|2009l ) to suggest a connection to 
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dusty environments in their host galaxies. We expect a 
correlation between variable absorption and the redden- 
ing within host galaxies, because formation of massive 
wind-swept shells at pc-scale distances from the SN can 
occur only in gas-rich hosts. But the line variability is 
confined to early times (Figure [3]), and the predicted de- 
cline is obviously in conflict with the much larger increase 
in N(Na I) observed in SNe2006X and 1999cl. 

Typical velocities of shells swept up by the progenitor 
WD winds are expected to be fairly low, ~ 10 km s~^, 
comparable to the intrinsic velocity dispersion of the dif- 
fuse ISM, because the SN explosion usually occurs in late 
stages of bubble evolution after the shell has been effec- 
tively decelerated by the ambient ISM. When the shell 
radial velocity falls within the saturated or partially sat- 
urated Na I D line profile of the host galaxy, prospects for 
the shell detection are markedly reduced. This strongly 
favors detection of blueshifted line absorption preferen- 
tially produced by progenitors with negative radial ve- 
locities (in the host galaxy frame), comprising no more 
than 50% of all SNe. Another important selection effect 
is associated with the ambient medium density; forma- 
tion of massive shells cannot occur in the low-density, hot 
ISM phase with a typical volume fraction of ~ 0.5. As a 
result, variable absorption is expected in only a fraction 
(certainly no more than a quarter) of all Type la SNe. 

Detection of variable Na I absorption from CSM shells 
is more difhcult than for SNe 2006X and 1999cl, and re- 
quires high-resolution spectroscopy soon after a SN dis- 
covery. But such observations can map the ambient 
medium around many Type la SNe, allowing for the de- 
tection of CSM shells and establishing their distances 
from the SN by measuring the rate of decrease in Na I 
absorption. Most progenitor scenarios predict the for- 
mation of shells at pc-scale distances in many progeni- 
tor systems. This prediction can be verified by obser- 
vations of variable absorption, advancing our knowledge 
of poorly understood Type la progenitors. We can learn 
much about these progenitors by correlating variable ab- 
sorption with SN properties such as luminosities, UV ex- 
cesses, or the presence or absence of high-velocity fea- 
tures. Observational efforts focused on high spectral res- 
olution observations of Type la SNe following soon after 
their discovery are warranted. 
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